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Interactions between non-BPS non-Abelian vortices are studied in non-Abelian 
U(l) X SU(N) extensions of the Abelian-Higgs model in four dimensions. In 
addition to the usual type I /II Abelian superconductors, we find other two new 
regimes: type I* /II*. 
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1. Introduction 

Recently, a new type of BPS vortex was found in U(N) gauge theo- 
ries. 1 ' 2 This is called non-Abelian vortex and carries the non-Abelian charge 

QpN-i _ S ^ N * >c J~tt/t\ ■ Readers can find good reviews in 3 ' 4 and 

references of related works therein. In this talk we are interested in studying 
interactions between non-Abelian vortices which are non-BPS. The non- 
BPS vortices are more natural than BPS ones in a sense that the BPS 
always requires a fine tuning or supersymmetry. It is well known that ANO 
vortices 5 ' 6 in the type I system feel an attractive force while those in the 
type II model feel a repulsive force. 7-10 Specifically we are interested in the 
interactions between vortices with different internal orientations, which is 
the distinct feature from the ANO case. 11 

This talk is based on 12 in collaboration with R.Auzzi and W.Vince. 



2. The model 

2.1. A fine-tuned model 

We start with non-Abelian, U (N) , extension of the Abelian-Higgs model 



C = Tr 



-^F^F^ + V^H^Hy - (v 2 l N - Hrff 



(1) 
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Here, for simplicity we take the same gauge coupling e for both the U(l) and 
SU(N) groups, while A 2 e 2 /4 is a scalar coupling and v (> 0) determines 
the Higgs VEV. H is N Higgs fields in the fundamental representation of 
U(N). The Higgs vacuum of the model is given by HW = v 2 1n. It breaks 
completely the gauge symmetry, although a global color-flavor locking sym- 
metry SU(N) c+f is preserved 

H - UqHUt, U g = Ul U G € SU(N) G , U F e SU{N) F . (2) 

The trace part TrH is a singlet under the color-flavor group and the traceless 
parts are in the adjoint representation. The U(l) and the SU(N) gauge 
vector bosons have both the same mass Mjj^ = Msu(N) = ev. The iV 2 
real scalar fields in H are eaten by the gauge bosons and the other N 2 (one 
singlet and the rest adjoint) have same masses M s = M a d = Xev. The 
critical coupling A = 1 (BPS) allows an TV = 2 supersymmetric extension. 



2.2. Models with general couplings 

A generalization of (1) is to consider different gauge couplings, e for the 
U(l) part and g for the SU(N) part, and a general quartic scalar potential 



Tr 



2^^ " h f ^ r +T>*H{V»H? 



V, (3) 



where we have defined = F*T A and /„„ = F° V T° with 

Tt{T a T b ) = S AB /2 and T° = 1/V2N The scalar potential is: 



\ 2 2 \ 2 2 

V = -^LTrX 2 + ^-Tr(X°T°-v 2 l N )\ (4) 

where = X°T° + X and X = 2£^J L 1 (H^T A Hi) T A . The symme- 

tries is same as the previous fine-tuned model (1). In this model, the U(l) 
and the SU(N) vector bosons have different masses = ev, M S u(N) — 

gv. Moreover, the singlet part of H has a mass M s different from that of 
the adjoint part M a( j as M s — X e ev, M a d = X g gv. For the critical values 
A e = \ g = 1 , the Lagrangian again allows an TV = 2 susy extension. 

2.3. Vortex equations in the fine-tuned model 

Let us make the following rescaling of fields and coordinates: 

H^vH, Wp^evWp, x^^. (5) 

ev 



2 
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The masses of vector and scalar bosons are rescaled to 

M u{1) = M SU(N) = 1, M s = M ad = A. (6) 

In order to construct non-BPS non-Abelian vortex solutions, we have 
to solve the equation of motion derived from the Lagrangian (1), 

VpF*™ - % - [H{V u Hy - {V U H)H^} = 0, (7) 

\2 

V^H + — (1 - HH ] ) H = 0. (8) 

From now on, we restrict ourselves to static configurations depending 
only on the coordinates x 1 ,^ 2 . Here we introduce a complex notation 
z = x 1 + ix 2 , = W = m ~j w * , V = 2i^2 = d + iW. In- 

stead of the equation of motions itself, it might be better to study gauge 
invariant quantities. For that purpose let us define 

W(z, z) = -iS^iz, z)BS(z, z), H(z, z) = S" 1 ^, z)H(z, z), (9) 

where S takes values in GL(N, C) and it is in the fundamental represen- 
tation of U(N) while the gauge singlet H is an N x N complex matrix. 
There is an equivalence relation (S, H) ~ (V(z)S, V(z)H), where V(z) is a 
holomorphic GL(N,C) matrix with respect to z. The gauge group U(N) 
and the flavor symmetry act as follows 

S(z, z) - U G S(z, z), H (z) - H (z)U F . (10) 

An important gauge invariant quantity is now constructed as 

Sl(z,z) = S(z,z)S(z,z)K (11) 

With respect to the gauge invariant objects, the equations of motion are 

4<9 2 (nan- 1 ) - ho (fftrr 1 ) + Bhh^- 1 = o, (12) 
no (n^Bn^j + B (pa (n- 1 ^ ) + — - hh^ n^H = 0. (13) 

These equations must be solved with the boundary conditions for k vortices 
detH — » z k , fl — » HH^ as z — » oo. 

2.4. BPS Limit 

For the later convenience, let us see the BPS limit A — > 1. It can be done 
by just taking a holomorphic function H with respect to z as 

H = H (z). (14) 
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Then the equations (12) and (13) reduce to the single matrix equation 

5 (non- 1 ) + 1 (l - HoHfo- 1 ^ = 0. (15) 

This is the master equation for the BPS non-Abelian vortex and the holo- 
morphic matrix Hq(z) is called the moduli matrix. 4,13 All the complex 
parameters contained in the moduli matrix are moduli of the BPS vortices. 
For example, the position of the vortices can be read from the moduli ma- 
trix as zeros of its determinant det H (zi) = 0. Furthermore, the number of 
vortices (the units of magnetic flux of the configuration) corresponds to the 
degree of det H (z) as a polynomial with respect to z. The classification of 
the moduli matrix for the BPS vortices is given in Ref. 4,13 

Consider U (2) gauge theory. The minimal vortex is generated by 

o-(l,0) _ (Z-Z 0\ tt(0-,1)_( 1 - b \ nf n 

H ° -{ -b> l)' H ° [oz-zoj- (16) 

zq corresponds to the position of the vortex and b and b' are the internal 
orientation. One can extract the orientation as the null eigenvector of Hq(z) 
at the vortex position z = z as 

* (i,o) =(i) - ^ (o ' i)= (i)- w 

Here "~" stands for an identification up to complex non zero factors: 
<j) ~ \(j>, A e C*, so that we have found CP 1 . 4,13 We call two non-Abelian 
vortices with equal orientational vectors parallel, while orthogonal orienta- 
tional vectors anti-parallel. 

Arbitrary two vortices (the center of mass is fixed to be zero and the 
overall orientaion is fixed) is given by 

The orientational vectors are then of the form 

3. Vortex interaction in the fine-tuned model 
3.1. (fci,fc 2 ) coincident vortices 

The minimal winding solution in the non-Abelian gauge theory is a mere 
embedding of the ANO solution into the non-Abelian theory. Embedding 
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is also useful for another simple non-BPS configurations. Let us start with 
the moduli matrix for a configuration of k coincident vortices. The axial 
symmetry allows a reasonable ansatz for and H 

n™" (J.?,))' " M = 0/(.W' (20) 

We call this "(0, A:)-vortex". When k > 2, it is possible that the ansatz (20) 
does not give the true solution (minimum of the energy) of the equations 
of motion (12) and (13). This is because there could be repulsive forces 
between the vortices. With ansatz (20) we fix the positions of all the vortices 
at the origin by hand. The master equation (15) is nevertheless still useful 
to investigate the interactions between two vortices. The results are listed 
in Table 1. For A = 1, the masses are identical to integer values, up to 10~ 5 



A 


k = 1 


k = 2 




0.8 


0.91231 


1.77407 


(o,a)-«oi»* 


0.9 


0.95737 


1.88936 


(.) typi I (A < 1) 


1 


1.00000 


2.00000 




1.1 


1.04053 


2.10655 




1.2 


1.07922 


2.20944 





Fig. 1. Spectrum of the (0, 2) and (1, 1) coincident vortices. 



order, which are nothing but the winding number of the vortices. 

There is another type of composite configuration which can easily be 
analyzed numerically 

(.i,i) = (Mr) ^ 0(i,i) = f/i(^ 1 

I Mr))' V f 2 {r)z^ 

This ansatz corresponds to a configuration with k\ composite vortices which 
wind in the first diagonal U(l) subgroup of U(2) and with ki coincident 
vortices that wind the second diagonal U(l) subgroup. We refer to these 
as a "(fei, fc2)-vortex". The mass of a (fci, fc2)-vortex is thus the sum of the 
mass of the (k\, 0)-vortex and that of the (0, fc2)-vortex. 

We call the non-Abelian vortices in the fine-tuned model for A < 1 type 
I, while they will be called type II for A > 1. From Fig. 1, we can see 
that in the type I case, the (0, 2)-vortex is energetically preferred to the 
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(1, l)-vortex, while in type II case the (1, l)-vortex is preferred. If the two 
vortices are separated sufficiently, regardless of their orientations, the mass 
of two well separated vortices is twice that of the single vortex. This mass 
is equal to the mass of the (1, l)-vortex. 

3.2. Effective potential for coincident vortices 

The dynamics of BPS solitons can be investigated by the so-called moduli 
approximation. 14 The effective action is a massless non- linear sigma model 
whose target space is the moduli space. If the coupling constant A is close 
to the BPS limit A = 1, we can still use the moduli approximation, to inves- 
tigate dynamics of the non-BPS non-Abelian vortices by adding a potential 
of order |1 — A 2 | <C 1. To this end, we write the Lagrangian 

C = £ B ps + (A2 ~ 1} (lAT-fffft) 2 . (22) 

We get non-BPS corrections of order 0(X 2 — 1) by putting BPS solutions 
into Eq. (22). The energy functional thus takes the following form 

£ = 2 + (A 2 -l)V, V = ^Jdx^x 2 Tr(l- |ff B ps(^)| 2 ) ' (23) 

where Hspsifi) stands for the BPS solution. We have defined a reduced 
effective potential V which is independent of A. The first term corresponds 
to the mass of two BPS vortices and the second term is the deviation from 
the BPS solutions which is nothing but the effective potential we want. 

To have the effective potential on the moduli space of coincident vortices, 
it suffices to consider only the matrix (18) with turning off the relative 
distance zq. In order to evaluate it, we need to solve the BPS equations 
with an intermediate value of rj. Because of the axial symmetry and the 
boundary condition at infinity il — ► Hq(z)Hq(z), we can make an ansatz 

The advantage of the moduli matrix formalism is that only three functions 
Wi(r) are needed and the formalism itself is gauge invariant. The effective 
potential can be obtained by plugging numerical solutions into Eq. (23). 
The result is shown in Fig. 2. 

The type II effective potential has the same qualitative behavior as 
showed in the figure. It has a minimum at \rj\ = 0. This matches the previous 
result that the (1, l)-vortex is energetically preferred to the (2, 0)-vortex. 
The type I effective potential can be obtained just by flipping the overall 
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Fig. 2. Numerical plots of the effective reduced potential V(|r?|). 



sign of that of the type II case. Then the effective potential always takes 
a negative value, which is consistent with the fact that the masses of the 
type I vortices are less than that of the BPS vortices. Contrary to the 
type II case, the type I potential has a minimum at \rj\ = oo, so that the 
(2,0)-vortex is preferred to the (1,1) vortex. 



3.3. Interaction at generic vortex separation 

Next we go on investigating the interactions of non-Abclian vortices in the 
U(2) gauge group at generic distances. We will again use the moduli space 
approximation. The generic configurations are described by the moduli ma- 
trices in Eq. (18). By putting the two vortices on the real axis, we can reduce 
zq to a real parameter d. So 2d is the relative distance and r\ the relative 
orientation. Now let us study the effective potential as function of r\ and d. 
As before, we first need the numerical solution to the BPS master equation. 
Despite the great complexity by broken axial symmetry, the moduli matrix 
formalism is a powerful tool and the relaxation method is very effective to 
solve the problem. Once we get the numerical solution, the effective poten- 
tial is obtained by plugging them into Eq. (23), see Fig. 3. It for the type 
II has the same shape, up to a small positive factor (A 2 — 1). The potential 
forms a hill whose top is at (d, \rj\) = (0,oo). It clearly shows that two 
vortices feel repulsive forces, in both the real and internal space, for every 
distance and relative orientation. The minima of the potential has a flat di- 
rection along the d-axis where the orientations are anti-parallel (77 = 0) and 
along the 77 axis at infinite distance (d = 00). Therefore the anti-parallel 
vortices do not interact. In the type I case (A < 1) the effective potential 
is upside-down of that of the type II case. There is unique minimum of the 
potential at (d, \n\) = (0, 00). This means that attractive force works not 
only for the distance in real space but also among the internal orientations. 
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Fig. 3. Left panel is the effective potential V(r],d). The Abelian potential V e 
(middle) and the non-Abelian potential Vg (right) for 7=1. 



4. Vortices with generic couplings 

In this section we sutudy the general model denned in Eqs. (3) and (4). We 
have three effective couplings 7 = g/e,X e ,X g after the rescaling (5). The 
masses of particles are rescaled as 



M u(i) = 1> Msu(N) = 1, M s = A e , M ad = 7A g 



(25) 



In order to find the effective potential on the moduli space as before, we 
need to clarify BPS configurations. The moduli matrix in (14) is still valid, 
while the master equation (15) get a modification 



49 (ftsrr 1 ) = ftofr 1 - i N + ( 7 2 - 1) ftofr 1 ^ — J -i N 



(26) 



where f2 = SS^ is same as before and Sl = H Hq. It turns out that the 
effective potential consists of the Abelian and the non-Abelian potentials 

V e (77,d; 7 ) = / dx 2 Tr(Fi 2 T ) 2 , V g ( V ,d;j) = [ di 2 Tr(F 12 ) 2 . (27) 



The true potential is a linear combination of them 



V(r],dn,Xe,X g ) = (Xl ~ l)V e (v,d;f) + 



Xi-1 



■V fl (»j,d;7). (28) 



4.1. Equal gauge coupling 7 = 1 revisited 



g — ^ e in the left panel of 



The effective potential with 7 = 1 and A = X g = A, 
Fig. 3 should be now decomposed in the two potentials, see the middle and 
the right panels in Fig. 3. In the case with X 2 — 1 > and X 2 — 1 > 0, the 
effective potential will have the same qualitative behaviors like the reduced 
potentials in the Figs. 3. The figures shows how V e and V g behaves very 
differently. In particular, the Abelian potential is always repulsive, both in 
the real and internal space. The non-Abelian potential is on the contrary 
sensitive on the orientations. Fig. 3 shows that it is repulsive for parallel 
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vortices while it is attractive for anti-parallel ones. When the two scalar 
couplings are equal, \\ = )? g , the two potentials exactly cancel for anti- 
parallel vortices. 

Of course, the true effective potential depends on A e and X g through 
the combination in Eq. (28). This indicates the interaction between non- 
Abelian vortices is quite rich in comparison with that of the ANO vortices. 



4.2. Different gauge coupling 7^1 

We now consider interactions between non-Abelian vortices with different 
gauge coupling e 7^ g (7 7^ 1). In Figs. 4 and 5 we show two numerical 
examples for the reduced effective potentials V e , V g given in Eq. (27). These 



anti-parallel (77 = 0) intermediate (77 = 4) parallel (77 = 00) 

Fig. 4. Effective potential with 7 = 1/2 vs. separation, (red, blue) = (Ve, V g ). 



anti-parallel (77 = 0) intermediate (rj = 4) parallel (77 = 00) 

Fig. 5. Effective potential with 7 = 1.3 vs. separation, (red, blue) = (V e , Vg) 



show that the qualitative features of V e and V g are basically the same as 
what is discussed in the equal gauge coupling case (7 = 1). The true effective 
potential in Eq. (28) depends on the three parameters 7, A e and X g . We can 
have potentials which develop a global minimum at some finite non zero 
distance, see Fig. 6 The figure shows the presence of a minimum around 
d ~ 2. This kind of behavior have not been found for the ANO type I/II 
vortices and the possibility of bounded vortices really results from the non- 
Abelian symmetry. 
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5. Interaction at large vortex separation 

5.1. Vortices in fine-tuned models e — g and A e = X g 

We study an asymptotic forces between vortices at large separation, follow- 
ing Refs. 10 We need to find asymptotic behaviors around (l,0)-vortex 

\{i,o)^l' z0 \ ,7(1,0) fl 



Ho{z) = {oi)> & [oj- (29) 

We are lead to the well known asymptotic behavior of the ANO vortex 

#[i.i] = (i + ^^o(Ar)) e\ = -i - £^i w) e\ (30) 

where K\ = —K' Q and we have defined -fffi i] and W^i] as [1,1] elements 
of H and W in Eq. (9) with the k = 1 ansatz (20). 

Next we treat the vortices as point particles in a linear field theory 
coupled with a scalar source p and a vector current j M . To linearize the 
Yang-Mills- Higgs Lagrangian, we choose a gauge such that the Higgs fields 
is given by hermitian matrix H = 12+5/jVj, = |iu* Cj with a = (I2, <?). 
with all h a , m)° are real. Then the quadratic part of the Lagrangian is 



•^free ^ ; 

a=0 



2 



(31) 



with /" = d^wf, — <3„io". We also take into account the external source 
terms to realize the point vortex 

3 

Source = [/^-iXl- ( 32 ) 
a=0 

The scalar and the vector sources should be determined so that the asymp- 
totic behavior of the fields in Eq. (30) are replicated. The solution of the 
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equation of motion is 



h° = h 3 = ±K (\r), p° =p z =qS(r), 



w°=w 3 = -^kxVX (r), j°=j 3 

ZTT 



-mk x V(5(r) (33) 



where k is a spatial fictitious unit vector along the vortex world- volume. The 
vortex configuration with general orientation is also treated easily, since the 
origin of the orientation is the Nambu-Goldstone mode associated with the 
broken SU(2) color-flavor symmetry H -> H (z)^U F , fa = U^' 0) . 
The interaction between a vortex at x = xi with the orientation fa and 
another vortex at x = X2 with the orientation fa is given through the source 
term and is summarized as 

2 



V w 



A fa 





2 




$1 




fa 



£*„(Ar) - £* (r) 



(34) 



where r = |xi — X2I 3> 1. When two vortices have parallel orientations, 
this potential becomes that of two ANO vortices. 10 On the other hand, 
the potential vanishes when their orientations are anti-parallel. This agrees 
with the numerical result found in the previous sections. In the BPS limit 
A = 1 (q = m), the interaction becomes precisely zero. 



5.2. Vortices with general couplings 

It is quite straightforward to generalize the results of the previous section 
to the case of generic couplings. We find the total potential Vi nt 
1 
2 



-K Q (X e r) H K (r) 



( 



+ 



2ir 




At large distance, the interactions between vortices are dominated by the 
particles with the lowest mass M\ ow . There are four possible regimes Vi nt = 




for M lo 



e -X 9 ir for MlQw = Mi 



ad • 



Type I 
Type I* 



for Afiow = Mu(i), Type II 
£- r e-rr for Af low = M su{2) , Type II* 



(36) 
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because of K (Xr) ~ yjn /2\re~ Xr . This generalizes the type I/II classi- 
fication of Abelian superconductors. We have found two new categories, 
called type I* and type II* , in which the force can be attractive or repulsive 
depending on the relative orientation. In the type I* case the forces be- 
tween parallel vortices are attractive while anti-parallel vortices repel each 
other. The type II* vortices feel opposite forces to the type I*. The result 
in Eq. (36) is easily extended to the general case of U(l) x SU(N). This 
can be done by just thinking of the orientation vectors <p as taking values 
in C^- 1 . 

It may be interesting to compare these results with the recently studied 
asymptotic interactions between non-BPS non- Abelian global vortices. 15 
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